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ABSTRACT 

The variability of the Fe Ka line near 6.5 keV seems to be reduced compared to 
the variability of the hard X-rays which presumably drive the line emission. This is 
observed both in active galactic nuclei and galactic black hole binaries. We point out 
that such reduced variability, as well as lack of coherence between the variations of the 
line and the continuum, are a natural prediction of a propagation model of variability 
in the geometry of inner hot accretion flow. We compute detail model predictions of 
the variability characteristics which could be compared with current and future data. 
We also point out that the model requires a gradual disappearance of the cold disc, 
rather than a sharp transition from the cold disc to a hot flow. 

Key words: accretion, accretion disc - binaries: general - X-rays: binaries - X-rays: 
galaxies - galaxies: active 



1 INTRODUCTION 

The Fe Ka fluorescent /recombination line near 6.5 keV is 
an important diagnostic of accretion flows around compact 
objects (see Reynolds & Nowak 2003 for a recent review). It 
is the strongest hard X-ray (E > 1 keV) line which can orig- 
inate in the innermost regions of accretion flows (< 100 R s ). 
It is indeed observed in energy spectra from all kinds of ac- 
creting sources: black hole and neutron star X-ray binaries, 
cataclysmic variables, active galactic nuclei (AGN). 

The line is produced when plasma is irradiated by hard 
X-ray (E > 7 keV) radiation. If the plasma is Thomson 
thick the continuum spectral component formed as a result 
of the irradiation ( "Compton reflection" ) has a characteris- 
tic shape, peaking at 20-30 keV (Lightman & White 1988). 
The line and the Fe K-shell absorption edge are superposed 
on this continuum (George & Fabian 1991; Matt, Perola 
& Piro 1991). The properties of the line and edge depend 
on the ionization of the reflecting medium: with increasing 
ionization the line and edge shift towards higher energies, 
while their strength increase (e.g. Zycki & Czerny 1994 and 
references therein). The profile of the line may be modified 
by relativistic and kinematic effects, if the line originates in 
e.g. a rotating accretion disc (Fabian et al. 1989). The line 
and absorption edge are then broadened and smeared. Such 
broad features are commonly seen in Seyfert galaxies (e.g. 
MCG-6-30-15, Tanaka et al. 1995, Fabian et al. 2002) and 
black hole binaries (BHB; e.g. Cyg X-l, Done & Zycki 1999; 
GRS 1915+105, Martocchia et al. 2002), providing a clear 
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evidence for a relativistic accretion disk extending deep into 
the gravitational potential of the central black hole. 

X-ray emission from accreting sources is highly vari- 
able in a broad range of time-scales. Most of the variability 
power is located in the range of Fourier frequency / w(0.1- 
1)M/(1OM ) Hz, corresponding to a time-scale (T = 1//) 
of a few seconds for a 10 Mq stellar black hole and a few 
weeks for a 10 7 M AGN (BHB review in McClintock & 
Remillard 2003; Markowitz et al. 2003b for AGN). Typical 
power density spectrum is roughly a power law with slope 
a ~ (P(f) oc f a ) at low /, steepening to a m —1 at 
/ w 0.1 Hz and to a « -2 at / = 1-3 Hz, for the BHB 
in low/hard state. Typical root-mean-square (r.m.s.) vari- 
ability is 20-30%. The variability is stochastic rather than 
caused by deterministic chaos type of process (Czerny & 
Lehto 1997). 

The observed variability of the Fe Ka line and the entire 
reprocessed component is somewhat surprising: they gener- 
ally show rather less variability than the high energy con- 
tinuum which is presumably driving the line emission and 
Compton reflection. This is seen both in AGN and BHB. In 
AGN the Ka line seems either not to respond to continuum 
variations on time-scales of minutes to days (e.g. Reynolds 
2000; Done, Madejski & Zycki 2000; Chiang et al. 2000), 
or the line variability appears to be uncorrelated with that 
of the continuum (Vaughan & Edelson 2001). In particular, 
studies of r.m.s. variability amplitude as a function of energy 
demonstrate the reduced variability in a relatively model in- 
dependent way (Inoue & Matsumoto 2001; Markowitz, Edel- 
son & Vaughan 2003a) . The short term variability of the re- 
processed component is more difficult to measure because of 
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insufficient statistics, but where such studies were possible, 
the results also suggested reduced variability (e.g. Done et 
al. 2000). 

A similar effect is seen in BHB systems. The time re- 
solved spectral analysis is more difficult for BHB than for 
AGN, but Fourier spectroscopy (study of Fourier-frequency 
resolved spectra; Revnivtsev, Gilfanov & Churazov 1999, 
2001) can be used to investigate the variability character- 
istics on short time scales. Indeed, the amplitude of the re- 
processed component in high Fourier frequency spectra is 
smaller compared to that in low Fourier frequency spectra 
(Revnivtsev et al. 1999, 2001). In other words, the repro- 
cessed component responds to the variability of the primary 
continuum on long time scales (T = 1// > 1 sec) but it does 
not do so on short time scales. 

In AGN one reason for the lack of variability might be 
the contribution to the reprocessing by a distant matter, e.g. 
the obscuring dusty torus. Indeed, a narrow Fe Ka line is 
observed in many Seyfert galaxies. However, it is the vari- 
ability properties of the broad component (hence produced 
very close to the central black hole) that are so puzzling. 

One suggested explanation invoked complex ionization 
effects in the illuminated surface of the disc (e.g. Nayak- 
shin, Kazanas & Kallman 2000). Formation of a hot ionized 
skin (where the Fe Ka line is not produced) with thickness 
proportional to the X-ray flux may lead to anti-correlation 
between the line flux and the continuum flux. This model 
was quantitatively tested by Zycki & Rozariska (2001), who 
concluded that it does indeed predict certain decrease of am- 
plitude of variability of the line, although it is not possible 
to obtain an absolutely constant line flux. 

In this paper we demonstrate that the reduced variabil- 
ity of the Fe Ka line is a natural result from a propagation 
model of X-ray emission in the geometry of a hot inner ac- 
cretion flow. This geometry is one of the possibilities for 
accretion flow in low/hard states of accreting black holes 
(e.g. Di Salvo et al. 2001; see Done 2002 for review). It is 
supported by X-ray spectral studies (review in Poutanen 
1999), in particular the correlation between spectral slope 
and amplitude of reflection (Zdziarski, Lubihski & Smith 
1999) and long time scale spectral evolution of black hole 
binaries (Zdziarski et al. 2003). The physical mechanism 
of formation the two-phase plasma flow may be related to 
plasma evaporation/condensation (Rozanska & Czerny 2000 
and references therein) . The idea of propagating X-ray emit- 
ting structures was put forward by Miyamoto et al. (1988), 
based on early variability studies of Cyg X-l. This was fur- 
ther developed and tested by e.g. Nowak et al. (1999), Misra 
(2000), and formulated in a more general form by Kotov, 
Churazov & Gilfanov (2001). Recently, Zycki (2003) showed 
that the Fourier-frequency resolved spectra can be repro- 
duced in the propagation model, while Uttley (2004) argued 
for this model based on the r.m.s.-flux relation in accreting 
pulsar SAX J1808.4-3658 (see also Uttley & M c Hardy 2001). 



2 THE MODEL 

The model was also described in detail in Zycki (2003). 
X-rays are assumed to be produced by compact active re- 
gions/emitting structures traveling from outside inwards, to- 
wards the black hole. The structures originate at a certain 



radius, R ou t, and they move towards the centre at a fraction 
of the free-fall speed, 

a a 1 2GM ... 
v = f3v B = /3d — — , (1) 

where /3 < 1. The plasma heating rate is assumed to depend 
on radius 

i h (t) oc R(t)- 2 b[R(t)], (2) 

where the exponent —2 corresponds to gravitational energy 
dissipation per ring of matter. Here I is the compactness pa- 
rameter, I = (L I ^flarc)o"T / (m c c 3 ) (where Dfl aro is the char- 
acteristic radius of the structure, assumed constant in time), 
and b(R) is the boundary term which is assumed to have its 
standard form b(R) = 1 — y^GRg/R (Shakura & Sunyaev 
1973). 

The heating rate is normalized to have an assumed max- 
imum value, Zh,max, the same for all flares. The equation of 
motion (Eq.0 can be solved to give 

r{t) = {rl£-Atr\ (3) 

where the radial positions r(t) and r out are expressed in units 
of R g = GM/c 2 (lowercase r will denote radial position in 
units of R g ). The duration of a flare is 

ttrav = (^ t 2 - rl( 2 )/A, (4) 

where the final radial position is assumed R ln = 6 R s . In ac- 
tual computations we treat i tr av as a parameter, and solve 
for (3. Following Poutanen & Fabian (1999) we generate ttrav 
according to a probability distribution P(t) oc t p for t be- 
tween i min and f max . 

Soft photons for Comptonization are assumed to come 
from reprocessing and thermalization of the hard X-rays, so 
that the feedback loop is realized as needed to explain the 
correlation between the amplitude of reflection and spectral 
index (Zdziarski et al. 1999; Gilfanov, Churazov & Revnivt- 
sev 2000). The geometrical scenario considered here is that 
of an inner hot flow partially overlapping a cold, optically 
thick disc disrupted at a certain radius, R tI , (see e.g. Pouta- 
nen, Krolik & Ryde, 1997, for arguments for the overlap). 
The luminosity of soft photons crossing an active region is 
parametrized as 

{1 for v Tt 

(Ji.y for r<Z, (5) 

where C(r) represents a covering factor of the cold repro- 
cessing matter. The covering factor is assumed to increase 
with R (i.e. 7 > 0), up to C(r tI ) = 1. The normalization 
constant iV s « 0.5, as appropriate for a continuous corona 
at r > rtt- Thus, during a flare the Comptonized spectrum 
has a constant slope of V « 2 up to the moment of crossing 
the truncation radius, and then the spectrum get harder (r 
decreases). The truncation radius and exponent 7 determine 
the average slope of the Comptonized spectrum. 

The emitting structures are assumed to originate at a 
mean rate of A per second. Time intervals between their 
launch are generated from the Aexp(— Xt) distribution, as 
appropriate for a Poissonian process. We adopt here the flare 
avalanches description (Poutanen & Fabian 1999), where 
each spontaneous ( "parent" ) flare has certain probability to 
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stimulate a number of "baby" flares. The stimulated flares 
are delayed after their parent flare (Poutanen & Fabian 
1999). 

Each emitting structure is followed until it reaches the 
marginally stable orbit at _R ms = 6 R g . At each time step 
the primary Comptonized component and its reprocessed 
component are computed. 

The Comptonized spectrum is computed using the code 
thComp (Zdziarski, Johnson & Magdziarz 1996), solving 
the Kompaneets equation. Computations are parametrized 
by the photon spectral index, which we compute from 
T = 2.33(ZhAs)~ 1/10 (Beloborodov 1999a,b), and electron 
temperature kT e /(m e c 2 ) computed using formulae from Be- 
loborodov (1999b). Plasma optical depth is assumed r cs = 1. 

The spectrum of the Compton reflected continuum is 
computed from the simple formula of Lightman & White 
(1988), 



S le g_(E) 



1 



~*Sprim(-£/) ; 



^abs ~\~ ^c: 



(6) 



where K, a b s (E) and K es (E) are the photo-absorption and 
electron scattering opacities, respectively, the former for a 
"cold" matter. The above formula is multiplied by a simple 
exponential cutoff to mimic the Klein-Nishina cutoff. The Fe 
Ka line is added to the reflected continuum, with equivalent 
width (EW) as a function of F, following computations of 
Zycki & Czerny (1994). Specifically, we find that the formula 
-0.5-0.15r 



EW(T) = 1.78 (r/1.1) 



keV 



(7) 



reproduces the dependence of EW (relative to the reflected 
continuum) on the spectral slope, F, for cold matter with 
solar iron abundance, and disk inclination of 30°. 

The instantaneous amplitude of the reprocessed compo- 
nent, fi/2-7r, should be related to the covering factor of the 
cold plasma, C(r) (Eq. |SJ. In practice we assume 



n/27r(r) = C(r) = 



for r > rtr 



fe) 7 for r < rtr 



(8) 



i.e. Q/2n decreases from 1 at r > r tr to Q./2n <C 1 at r = r ms . 

The sequence of spectra created by the above procedure 
is subject to standard analysis in the time and Fourier do- 
mains (see e.g. van der Klis 1995; Nowak et al. 1999; Pouta- 
nen 2001). 

All computations are performed assuming the central 
black hole mass M = 1O 7 M0. The time-scales of variabil- 
ity are assumed to follow a simple scaling with the black 
hole mass. Consequently, the values of parameters of the 
variability model are adjusted so that the power spectrum 
density (PDS) has the same shape as PDS of Cyg X-l in 
low/hard state, and is shifted in frequency by the mass ra- 
tio. This is achieved for i m i n = 5 x 10" 3 Mi SeC, tmax — 2iWl 
sec, p = —1, where Mi = M/(10Mq). The emitting regions 
are launched at R ou t = 100 R g at a rate A = 40Mj~ sec - . 
The truncation radius is assumed R tI = 30i? g , while 7 = 2. 
We consider a stationary model, i.e. the above values are 
constant in time. 



3 RESULTS 

The flux in the Ka line (the line light curve) is found simply 
by subtracting the continuum contribution from the total 
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Figure 1. During a flare the Ka line flux peaks earlier than 
the 7-30 keV continuum flux. Since the line flux can be written 
^line(i) oc f2/27r(t) X F(t), it decreases when the amplitude of 
reflection, Q/2iz, decreases for r < rtr- Note that the line and 
continuum fluxes were rescaled to 1 at the maximum, so the flare 
amplitudes in these two energy bands are not represented accu- 
rately. 



flux in the line bin. The Fe Ka line predicted by our model 
is narrow because of relatively low velocities of the active 
regions. This is required to reproduce the relatively long 
time-scales of maximum power, but it is also a consequence 
of the assumption that the active regions move purely radi- 
ally. If they possessed a significant azimuthal velocity com- 
ponent, their trajectories would be longer spirals and the 
total velocity would have to be accordingly higher, in order 
to produce the same flare durations. The azimuthal velocity 
component would in that case give rise to Doppler broaden- 
ing of spectral features, but the variability properties of the 
model would remain unchanged. 

Temporal profile of a typical flare is plotted in Fig. 
The profile of the heating rate, lh(t), and total Comptonized 
emission are similar to the plotted profile of the 7-30 keV 
flux. The amplitude of reflection, fl/2-K, remains constant 
at 1 until the emission region has reached the truncation 
radius, r tr , and then Q/2n begins to decrease. This has an 
important consequence for the flux of the Ka line: the line 
flux peaks earlier than the continuum flux. This is simply 
because the line flux can be expressed as the product of 
Q/2tt and -Fcont, Fu nc oc £l/2n x F cont . The magnitude of 
the difference in peaks position depends on the exponent 7 
describing the Q/2n(r) dependence (Eq. [SJ, and it is ~ 0.2 
of the total flare duration for the assumed 7 = 2. The second 
important consequence of the obtained Sl/2n(t) relation is 
that the amplitude of the flare in the Ka line is smaller than 
in the continuum (this is not shown in Fig. . 

Both the above effects have important consequences for 
line and continuum light curves, examples of which are plot- 
ted in Fig. |21 Variations of the line are clearly smoother that 
those of the continuum: there is less high-frequency variabil- 
ity and the peak-to-peak amplitude is clearly smaller in the 
line light curve, although major events in the continuum are 
reflected in the line light curve. The decoupling of variabil- 
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Figure 2. Part of a typical light curve: 9 keV continuum and the Ka line. Panel (a): 10 5 sec (~ 1 day), panel (b): 10 6 sec. On time scale 
the variability of the line is clearly smoother (less high-frequency variations, smaller peak-to-peak amplitude) than of the continuum. 
This is still visible, although less so on the longer time scale plotted. The black hole mass assumed for scaling the time scales is 10 7 Mq. 



ity is more clearly visible on the short term light curve (10 5 
sec), the differences between the line and continuum curves 
decreasing for longer time-scales. The curves are very similar 
when viewed on time scales of order ~ 10 7 sec. Another way 
of demonstrating the reduction of variability is to compute 
the power spectrum. This is plotted in Fig. [3J together with 
a number of other characteristics. The PDS of the line cuts 
off above /i w 10~ 6 Hz much more rapidly than PDS of the 
continuum. The r.m.s. computed for / > fi is 0.28 for the 
continuum, while it is only 0.16 for the line. The frequency 
/i is related to the time-scale of longest individual flares 
(tmax = 2 x 10 6 sec in our computations). 

There is a clear lack of coherence between the line and 
continuum light curves, as shown in panel (b) of Fig. [3] The 
coherence plotted there is defined as in Vaughan & Nowak 
(1997), and it is a measure of linear relation between two 
light curves. Generally, observed light curves show coher- 
ence very close to 1 in a broad range of Fourier frequency 
(Nowak et al. 1999 for Cyg X-l; Vaughan, Fabian & Nan- 
dra, 2003, for MCG-6-30-15). This indicates that each active 
region produces radiation in the entire X-ray band, and so 
photons of different energies track each other's variations 
very closely. Here, however, the line photons do not follow 
the continuum photons in a simple way, because of the vary- 
ing geometrical factor C(r). Hence the transformation be- 
tween the continuum and line photons is complicated and 
the coherence is reduced below unity. 

The time delay between the line and continuum pho- 
tons is much longer than between the two continuum bands 
of similar energy separation (the line is leading the variabil- 
ity of the continuum at E > 7 keV). Figure^ demonstrates 
this effect. Time lags between the 6 keV and 9 keV contin- 
uum bands would be shorter than the plotted time lag be- 
tween 3 keV and 9 keV bands. However, the time lag between 
the 6.4 keV Ka line and the 9 keV band is actually longer 
(by a factor of « 5) than the latter. This is a direct conse- 
quence of the line flux peaking earlier than the continuum 
during a flare, as already discussed (Fig.0. For example, in 



a ~day-long observation (Fourier / = 10 -5 -10 -4 ), the line 
may lead the 9 keV continuum by as much as 1-3 hours. The 
cross-correlation function between the Ka and 9 keV light 
curves is asymmetric, peaking at At w — 10 3 sec, consistent 
with the time lags. The cross-correlation functions plotted, 
CCF(t), were normalized by dividing them by the square 
root of the product of variances, which means that for a per- 
fectly correlated signals CCF(0) = 1. We can then observe 
that while two continuum bands are indeed well correlated, 
the correlation between line and continuum light curves is 
weaker, resulting in CCF(0) w 0.82. 



The model r.m.s. spectra show a minimum at the en- 
ergy of the Ka line, as indeed observed in spectra of Seyfert 
galaxies (Markowitz et al. 2003a) . The minimum is of course 
a direct consequence of the weakly variable line flux, which 
contributes to the count rate in that bin, but only weakly 
so to the r.m.s. variability. Since in our model the Ka line 
is narrow, the reduction of r.m.s. appears only at the line 
energy bin. Observed lines are broad, therefore the observed 
r.m.s. (E) dependence have broad minima around 6-7 keV. 



Fig. shows the zero-lag correlation between the 7— 
30 keV flux and the Ka line flux. Each point represents a 
spectrum averaged over ~ 5000 sec, meant to roughly repre- 
sent one RXTE orbit. When the two fluxes are low they are 
simply linearly related. This corresponds to emission from 
the outer regions, r > r tr , where the amplitude of reflection 
is Q./2-n = 1. Higher fluxes correspond to emission coming 
at least partially from the inner region, below r tr , where 
Q./2-k < 1, and so the linear relation breaks down. Fig. Et> 
shows the F(7 — 30) vs. F(Ka) for flux-binned spectra from 
the same, 10 7 sec light curve. In this representations the 
spatial separation discussed above is averaged over to some 
extent, resulting in significantly weaker relation, with loga- 
rithmic slope of « 0.6. 
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Figure 3. Variability properties of the Ka line vs. those of the 9 keV continuum, (a) power spectra, (b) coherence function, (c) time 
lags and (d) cross-correlation function. The Ka line shows weaker high-frequency variability compared to the continuum and there is a 
marked loss of coherence between the line and the continuum, both effects occurring above the frequency corresponding to longest flares. 
The time lags of the 9 keV continuum relative to the line are much longer than those relative to the 3 keV photons (factor of ~ 5), 
despite smaller energy separation. The lags can also be seen from the highly asymmetrical cross-correlation function. 



4 DISCUSSION 

We have computed the variability properties of the Fe Ka 
line from X-ray reprocessing in a propagation model of X- 
ray emission in accreting compact objects. The model com- 
bines results from various spectral and timing studies of 
accreting black holes. The former suggest the geometry of 
the standard optically thick accretion disc truncated at a 
radius larger than the radius of the last stable orbit (e.g. 
Esin, McClintock & Narayan 1997; Gierlihski et al. 1997; 
Zycki, Done & Smith 1998; Done & Zycki 1999; Zdziarski 
et al. 1999, 2003). The latter postulate correlated flares 
(avalanches) and spectral evolution during flares, in order 
to explain power spectra and time lags (Poutanen & Fabian 
1999; Kotov et al. 2001). Combined spectral-timing stud- 
ies suggest a connection between the geometry and timing 
properties (Revnivtsev et al. 1999, 2001; Zycki 2002, 2003). 
Most of the observational results were obtained for black 
hole binaries, but, were possible to conduct, analysis of data 
for Seyfert galaxies confirm the general similarity between 
the two classes of objects (e.g. spectral studies by Done et 
al. 2000; Chiang et al. 2000; Lubiriski & Zdziarski 2001; tim- 



ing studies of Czerny et al. 2001; Uttley & M c Hardy 2001; 
Vaughan et al. 2003; Markowitz et al. 2003b; Papadakis 
2004). 

The variability properties of the Fe Ka line seem to be 
similar in both classes of objects in that the line appear to 
be less variable than the continuum which drives it, and, 
where the variability is detected, it does not seem to be 
clearly correlated with the continuum. This is contrary to 
simple(st) ideas, whereby the continuum and the line are 
produced in the same region and thus should be closely re- 
lated. We note that the observational situation is far from 
clear, though. Time resolved spectral analysis would be the 
most direct method to determine the variability of the line, 
but the results may be model dependent (see detailed discus- 
sion in Zdziarski et al. 2003) . Decomposition of counts from 
a medium energy resolution instrument like RXTE/PCA 
into the line and continuum depends on the model assumed 
for both the line and the reflected continuum. In particular, 
the effects of relativistic smearing of the reflected continuum 
was not taken into account. We note though that according 
to Markowitz et al. (2003a) the results on line variability in 
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Figure 4. Dependence of r.m.s. variability on energy. Clear dip 
at Fe Kq line energy is seen, from the weakly variable line. 



a sample of Seyfert galaxies were insensitive to assumptions 
about the amplitude of the reflected component: whether 
its relative amplitude was fixed, or allowed to vary in ac- 
cord with the Ka line. 

Our model does reproduce the reduced variability of the 
Fe Ka line, compared to the variability of its driving con- 
tinuum. Line variations may also appear not exactly corre- 
lated with continuum variations, because of the time delay 
between the peaks of the line and continuum fluxes. Both 
effects are necessary consequences of the adopted geometry, 
of a truncated disc with inner hot flow. The same geome- 
try can also explain the hard X-ray time lags (Kotov et al. 
2001; Zycki 2003) through spectral evolution during flares 
(Poutanen & Fabian 1999). Quantitatively, our assumed ra- 
tio of heating to cooling rates, C(r) oc lh(r)/l a (r), is a much 
weaker function of radius than could be expected for a com- 
pact (size <JC r) active region and a sharply truncated disc. 
As already mentioned in Zycki (2003), in that latter geom- 
etry the supply of soft photons from the disc would dimin- 
ish so rapidly that the predicted energy spectrum would be 
much too hard to be consistent with the data. This implies 
a gradual disappearance of the cold disc, which may be an 
interesting clue as to how the physical process of disc evap- 
oration proceeds (Rozariska & Czerny 2000). 

We emphasize that the reduction of variability of the 
Ka line discussed in the present paper is the same phe- 
nomenon as the decreasing reflection amplitude with Fourier 
frequency, found in X-ray data of Cyg X-l and GX 339-4 by 
Revnivtsev et al. (1999, 2001) and modelled by Zycki (2003). 
In the present paper this effect was analysed with tools usu- 
ally applied to AGN data, which are usually analysed in time 
domain rather than Fourier domain. We note, that it does 
not seem possible to design the parameters of the model, so 
that the line flux remains exactly constant. This is because, 
even though the line may be constant during each flare (as- 
suming C(r) oc 1/Zh(r)), the number of flares active at any 
time varies, and this causes variation of the total flux of the 
line. 

A robust feature of the presented model is a connection 
of the time scale of the line response to the time scale of 
continuum variability. This is because the line flux responds 
on the time scale related to the duration of a flare. This in 



turn has to be chosen such that the observed power spec- 
tra are reproduced. The peaks in / x P(f) are at rather 
long time-scales, much longer than just the light travel time 
through the region of most efficient energy generation. For 
example, the longest flares in our computations last w 10 6 
sec, which in light travel time corresponds to large distance 
of 2 w xl0 4 i? g (for 10 7 M Q ). This is simply a manifesta- 
tion of the well known observational fact that time-scales 
of maximum X-ray variability power are much longer than 
the naively expected short dynamical time-scale. The phys- 
ical mechanisms of producing the relatively long time scales 
of variability are rather unclear, but some interesting pos- 
sibilities were recently considered in literature. 3-D magne- 
tohydrodynamical simulations of Narayan, Igumenshchev & 
Abramowicz (2003) reveal slow (v <C vb) drift of plasma 
clumps across lines of magnetic field, which lines are com- 
pressed by the accretion flow. King et al. (2004) considered 
a model involving magnetic dynamos operating locally in 
the disc. The longer and larger flares are results of cor- 
relations between dynamos acting in neighboring locations 
(radii). While individual dynamos operate of short (dynam- 
ical) time scales producing short flares, the correlated be- 
haviour produces longer lasting, large events. Whatever the 
exact physical processes are, it is clear that the accretion 
flow is highly inhomogeneous and structured, factors that 
any realistic modelling should allow for. 

In the presented model, the line flux is leading the con- 
tinuum flux. We ignored the light travel time delay of the 
line photons after the continuum, but this is unlikely to af- 
fect our result, since the line is supposed to originate close 
to the location of emission of primary radiation. The light 
travel time delay may be estimated as St ~ r x R s /c (as- 
suming the height of the emission region h ~ r), which gives 
St ~ 10 3 sec, i.e. about the time bin in our simulations. Ob- 
viously, there may be additional effects due to, for example, 
adjustment of properties of the reprocessing medium to the 
increasing irradiation flux, which might affect the result to 
some extent. 

The dependence of r.m.s. variability amplitude on en- 
ergy is a relatively model independent demonstration of the 
reduced variability of the line. Our computations qualita- 
tively reproduce the minimum of r.m.s. (E) at the energy 
of the line. Generally, the r.m.s. spectra are energy depen- 
dent variability amplitude, or, equivalently, can be thought 
of as representing energy spectra of the variable component 
of the spectrum. This can in principle be computed also 
in narrow ranges in Fourier frequency (Fourier frequency 
resolved spectroscopy, Revnivtsev et al. 1999, 2001; Zycki 
2002, 2003), but the quality of AGN data is not sufficient to 
make use of this technique as yet. 

Concluding, the propagation model of high energy emis- 
sion in the geometry of a truncated accretion disc provides a 
framework for understanding many of the observed spectral 
and temporal characteristics of X-ray radiation from accret- 
ing black holes. 
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Figure 5. Line flux vs. the 7-30 keV flux. Left panel: each point shows one time bin of 5000 sec. from a 10 7 sec. observation. The 
line flux follows linearly the continuum flux for low values of the latter, corresponding to the emission coming from outer regions (active 
region above the cold disc). The higher the continuum flux, the stronger the deviation from a linear relation, and the bigger the spread 
of points. Right panel: results from spectra binned in flux. The slope of the best fit line is ?s 0.56. 
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